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摘 要 : 降水 量 作 为 人 工 增 水 效果 分 析 常 用 指标 , 受 地 理 、 经 济 、 技 术 等 影响 ,往往 可 研究 分 析 的 代 
表 性 数据 站 点 数量 较 少 ,一 定 程度 上 对 区 域 性 效果 检验 精度 造成 影响 。 为 此 ,基于 1973 一 2018 年 
的 5 一 9 月 巴 音 布鲁克 气象 站 逐日 气象 资料 和 开 都 河上 游 大 山口 水 文 站 逐 月 径流 量 资料 ,利用 Bu- 
dyko 模 型 构建 径流 模拟 方程 ,并 运用 序列 试验 检验 、 不 成 对 秩 和 检验 以 及 1 检验 等 统计 方法 ,以 降 
水 量 、 径 流量 等 作为 指标 ,探讨 该 区 域 暖 季 不 同 统计 指标 对 人 工 增 水 作业 效果 检验 的 差异 性 。 结 


果 表 明 :(1) 基于 Budyko 模 型 得 出 的 径流 量 与 降水 


量 不 但 相关 性 极 高 (及 =0.9971,P<0.001) ,而 且 


速率 与 趋势 变化 一 致 ,表明 模拟 径流 量 不 仅 能 准确 反映 降水 量变 化 趋势 ,还 能 代表 降水 对 径流 的 
影响 量 ,(2) 将 实测 径流 量 、 模 拟 径流 量 和 降水 量 作为 统计 变量 ,利用 不 对 称 秩 和 检验 及 t 检 验 , 分 
析 得 出 人 工 增 水 作业 后 ,降水 量 与 径流 量 增 加 显著 (P<0.02)。(3) 降水 量 作为 统计 指标 检验 功效 最 
好 ,人 工 增 水 作业 开展 后 只 需 增值 11.59% 就 能 显著 检验 出 效果 ,而 模拟 径流 量 相 比 实测 径流 量 的 检 
验 功 效 值 低 3.72% ,说 明 检验 功效 有 提升 。(4) 选取 统计 显著 性 水 平 90% 的 置信 区 间 , 得 出 作业 期 
(1994 一 2018 年 ) 比 历史 期 (1973 一 1993 年) 的 暖 季 月 均 降 水 量 绝 对 增加 值 为 5.38 mm, 相 对 增 率 为 
12.05%; 模 拟 径流 量 绝对 增加 值 为 4.53 ms ,相对 增 率 为 14.7%; 实 测 径 流量 绝对 增加 值 为 28.48 mts”, 
相对 增 率 为 18.48% ,表明 巴 音 布鲁克 山区 暖 季 人 工 增 水 作业 效果 显著 。 
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文章 编号 : 1000 - 6060(2023) 12 - 1963 - 10( 1963 ~ 1972) 


水 资源 是 人 类 社会 可 持续 发 展 的 重要 物质 保 
障 之 一 ,可 随 着 人 口 增长 .社会 发 展 ,人 类 对 水 资 
源 的 需求 越发 加 剧 。 根 据 美国 科学 家 Schaefer 和 
Vonnegut 相继 提出 向 冷 云 播撒 催化 剂 ( 干 冰 ta 
银 ) 可 促使 降水 形成 的 办 法 ” ,一 些 国家 纷纷 兴起 人 
工 增 水 试验 研究 ”, 这 为 人 工 影响 天 气 发 展 及 解决 
水 资源 短缺 方法 提供 了 重要 指导 作用 “”。 我 国 作 
为 世界 的 贫 水 国之 一 ,人均 水 资源 占有 量 不 足 世 
界 人 均 1/4, 水 资源 的 资源 豪 赋 度 在 时 空 上 更 存在 明 
显 差异 ” ,特别 是 西北 地 区 ,水 资源 已 是 限制 当地 经 
济 社会 实现 高 质量 发 展 的 重要 因素 "”。 因 此 ,在 我 
国 西北 干旱 地 区 利用 人 工 增 水 技术 不 断 开 发 空中 
云 水 资 源 的 活动 ,已 成 为 各 地 人 为 增加 空中 降水 的 
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重要 措施 之 一 ,并 在 当地 生态 环境 保护 与 修复 .水 
资源 安全 保障 等 方面 发 挥 着 重要 作用 。 

随 着 我 国 各 级 政府 和 社会 的 大 力 支持 ,人 工 增 
水 作业 的 需求 和 规模 不 断 增 长 ,而 增 水 作业 后 的 效 
果 情 况 越发 受到 各 界 关注 ,可 自然 云雨 变化 复杂 ， 
如 何 科 学 合理 评 佑 人工 增 水 作业 效果 依然 困难 。 
目前 ,只 有 以 色 列 在 20 世 纪 60 一 70 年 代 开 展 的 人 
工 增 雨 催化 试验 得 到 广泛 认可 ,其 增 雨 效果 为 6%~ 
15%"'"。 较 为 著名 的 试验 还 有 美国 怀俄明 人 工 影 响 
天 气 试验 项 目 (WWMPP), 取 得 每 年 10%~15% 的 增 
雪 作 业 效 果 ”。 而 国内 最 著名 的 福建 古田 水 库 人 
工 增 雨 试验 ,通过 科学 化 作业 , 增 雨 效果 达 20% 以 
上 3。 近 些 年 ,学 者 们 运用 多 种 方法 进一步 研究 了 
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人 工 增 水 作业 的 机 理 与 效果 ,如 : 刘 平 等 “基于 雨 
滴 谱 观测 数据 开展 人 工 增 十 作业 微 物理 特征 分 析 ， 
表明 人 工 增 十 作业 后 ,大 滴 粒 子 显 著 增 多 ,雨滴 谱 
变 宽 ;李斌 等 利用 区 域 回 归 分 析 法 对 火箭 增 水 作 
业 效 果 进 行 研究 ,得 出 地 面 火 箭 人 工 增 水 作业 效果 
与 作业 方法 间 的 关系 ; 陈 宝 君 等 “利用 数值 模拟 方 
法 讨论 了 强 对 流 云 实 施 人 工 催化 作业 后 ,降雨 量 得 
到 有 效 增加 的 情况 。 以 上 研究 表明 ,科学 开展 人 工 
增 水 作业 是 具有 一 定 成 效 的 ,而 运用 合适 的 分 析 方 
法 ,能 有 效 检验 出 增 水 效果 。 

新 疆 位 于 我 国 西北 干旱 区 , 因 “ 三 山 夹 两 分 "的 
地 理 影 响 , 形 成 了 水 资源 山区 多 ,平原 少 的 特点 "， 
为 缓解 水 资源 短缺 ,开发 山区 空中 云 水 资源 已 成 为 
当地 重要 措施 之 一 。 作 为 我 国 西北 干旱 区 最 大 降 
水 中 心 ,被 称 为 干旱 区 ”* 湿 岛 " 的 天 山 山 区 ,一 直 
是 新 疆 人 工 增 水 作业 重要 地 区 ,而 位 于 其 腹地 的 巴 
音 布 鲁 克 山区 ,不 但 是 新 疆 经 济 发 达 区 域 的 主要 水 
资源 供应 地 之 一 ,还 是 天 山 最 大 的 人 工 增 水 作业 
区 。 然 而 ,由 于 地 理 、 经 济 .技术 等 因素 , 巴 音 布 鲁 
克 山 区 气象 探测 手段 有 限 , 人 工 增 水 作业 前 后 云雨 
变化 无 法 有 效 监测 ,以 及 效果 检验 的 要 求 较 高 ”， 
该 区 域 的 人 工 增 水 作业 效果 检验 研究 还 鲜 有 。 根 
据 已 有 研究 结果 ”” ,选择 人 工 增 水 作业 区 域 相关 
的 径流 量 作为 统计 变量 ,也 能 科学 检验 出 增 水 效 
果 ,考虑 到 地 表 径 流 形成 因素 较为 复杂 ,直接 选 
择 实测 径流 量 分 析 ,会 对 检验 结果 产生 较 大 误 
差 。 为 此 ,本 文 基于 Budyko 水 热 耦合 平衡 假设 方 
程 模拟 出 受降 水 影响 的 长 时 序 径流 量 数据 ,以 期 
提高 效果 检验 精度 ,并 采用 适当 统计 学 方法 对 巴 


z 


音 布鲁克 山区 暧 季 多 年 地 面 人 工 增 水 作业 效果 进 
行 对 比分 析 。 


1 数据 与 方法 


1.1 研究 区 概况 

巴 音 布鲁克 山区 位 于 新 疆 天 山中 上段 南 芳 , 由 大 、 
小 尤 尔 都 斯 2 个 高 位 山 间 盆地 和 山区 丘陵 草场 组 
成 ,其 内 不 仅 开 都 河 支流 遍布 ,还 有 一 定 面积 的 永久 
性 积 雪 和 冰川 (图 1)。 该 区 域 海拔 高 度 普遍 在 2400 m 
以 上 ,高 寒气 候 显 著 , 全 年 气温 较 低 , 蔡 发 量 较 少 , 冬 
长 夏 短 , 积 雪 日 数 达 139.3 d; 唯 有 了 暧 季 5 一 9 月 ,是 全 
年 最 适宜 开展 人 为 活动 的 时 段 ,期 间 月 均 气 温 相 对 
稳定 在 10 上 下 ,并 呈 逐 年 上 升 趋势 的 ,同时 全 年 
80% 的 降水 量 与 降水 日 数 和 年 径流 量 主要 分 布 时 段 
也 集中 在 该 时 期 。 为 此 , 巴 音 布 鲁 克 山区 每 年 
5 一 9 月 大 力 开 展 人 工 增 水 作业 ,以 期 缓解 日 益 加 剧 
的 水 资源 短缺 问题 ,如 今 该 区 域 已 成 为 新 疆 巴 音 郭 
榜 蒙 古 自治 州 最 大 的 人 工 增 水 作业 区 。 
1.2 数据 来 源 

所 需 资料 包括 日 降水 量 \ 日 气温 (平均 /最 高 /最 
低 )、 相 对 湿度 、 日 照 时 数 和 风速 等 日 气象 数据 以 及 
实测 月 径流 量 数据 ,分别 来 源 于 巴 音 布鲁克 国家 基 
本 气象 站 和 开 都 河流 域 出 山口 控制 站 (大 山口 水 文 
站 ) ,时 间 序 列 为 1973 一 2018 年 5 一 9 月 。 根 据 新 疆 
巴 音 郭 楞 蒙 古 自治 州 人 工 影 响 天 气 办 公 室 的 人 工 
增 水 作业 数据 显示 , 巴 音 布鲁克 山区 人 工 增 水 作业 
自 1994 年 开始 ,每 年 5 一 9 月 为 人 工 增 水 作业 主要 
时 有 段 ,期 间 年 均 作 业 300 次 ,作业 量 500 枚 。 为 方便 
后 续 研 究 , 将 未 开展 人 工 增 水 作业 的 1973 一 1993 年 
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图 1 巴 音 布鲁克 山区 站 点 分 布 示意 图 


Fig. 1 Schematic diagram of the distribution of stations in Bayanbulak mountain area 
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命名 为 历史 期 ;开展 人 工 增 水 作业 的 1994 一 2018 年 
命名 为 作业 期 。 所 有 气象 水 文 数据 都 严格 经 过 人 
工分 析 审 核 ,确保 数据 质量 控制 。 
1.3 研究 方法 
1.3.1 基于 Budyko 假 设 解析 表达 式 的 径流 量 模拟 
方程 由 于 巴 音 布鲁克 山区 地 表 径 流 主 要 补给 以 降 
水 和 冰雪 融 水 为 主 ,为 避免 检验 结果 产生 较 大 误 
差 ,根据 能 定量 得 出 水 文 循 环 中 降水 量 和 潜在 蒸 散 
发 量 对 径流 影响 程度 的 Budyko 水 热 耦 合 平衡 假设 
方程 ,以 及 独立 流域 长 期 下 垫 面 未 发 生变 化 时 ， 
其 径流 量 可 视 作 降水 量 与 实际 莱 散 发 量 之 和 的 水 
量 平 衡 关 系 当 ,构建 出 研究 区 主要 受降 水 影响 的 径 
流 模拟 方程 。 具体 公 式 如 下 : 
R's (pr + ET Y -ET, (1) 

式 中 : R' 为 受降 水 量 影响 的 径流 量 (m…s"'); w N 
个 积分 常数 ;P 为 降水 量 (mm) ;ET 为 潜在 燕 散 量 
(mm*d"')。ET。 采 用 世界 粮农 组 织 推 荐 的 通用 公式 
FA056-PM 计 算得 出 *。 
1.3.2 序列 试验 法 序列 试验 检验 法 是 人 工 增 水 试 
验 中 常用 的 一 种 效果 检验 方法 ,主要 是 将 统计 变 
量 的 历史 期 平均 值 作为 作业 期 的 自然 期 待 值 ,再 与 
作业 期 实测 值 比较 ,确定 增 水 效果 。 为 此 ,通过 比 
较 21 a 历史 期 与 25 a 作业 期 的 各 统计 变量 平均 值 差 
值 变化 ,分 析 人 工 增 水 效果 。 
1.3.3 不 成 对 秩 和 检验 法 不 成 对 秩 和 检验 法 是 一 
种 非 参 数 性 检验 5 ,可 对 人 工 增 水 作业 前 后 统计 变 
量 的 变化 情况 进行 显著 性 检验 。 其 原理 为 :2 个 总 
RACIB C) ,分 别 独 立 选取 2 个 容量 为 n, 和 nn 的 
随机 样本 (np) :away t, buo 4 
A(x)=B(x) ,接受 原 假设 Ho:2 个 总 体 相等 ;否则 接受 
假设 Hi, 即 人 工 增 水 效果 显著 。 一 般 先 将 统计 变量 
按照 一 定 顺 序 进行 排列 计算 出 秩 和 (7) ,然后 选取 
样本 容量 小 的 7 做 比较 。 

在 计算 中 , 当 n 与 分 别 大 于 10 时 ,7T 近 似 于 正 


态 分 布 N [mite — a 1) AA $ D) HA n, 


为 选取 7 的 样本 容量 。 此 时 ,统计 检验 值 (w) 可 通过 
FREH: 
n(n, +n, + 1) 
yo Pa Salton 7 yy 
BÆ 


nn(n, Eyt 1) 
12 


当 双边 检验 时 心 值 落 在 (-1.96,+1.96) 之 内 , 2 
异 不 显著 ; 若 v 值 落 在 (-1.96,+1.96) 之 外 ,差异 显 
著 , 显 闭 性 水 平 为 0.05。 单 边 检 验 时 , 若 x1.64( 或 uv 和 
-1.64) , 则 差异 显著 ;否则 不 显著 ,显著 性 水 平 o=0.05。 
1.3.4 1 检验 法 当 统 计 变量 在 总 体 服从 正 态 分 布 ， 
此 时 检验 分 布 参数 ,可 采用 参数 性 检验 。 当 样本 数 
量 z<30 时 ,参数 性 检验 选用 :检验 法 最 适当 ,但 要 求 
统计 变量 的 方差 不 改变 小 。 因 此 ,选择 适用 小 样本 
数据 正 态 检验 的 Shapiro- Wilk (S-W) 检 验 法 和 方差 
检验 的 下 检验 法 ,对 样本 的 正太 分 布 和 方差 齐 性 情 
况 进 行 显著 性 检验 。 其 中 ,S-W 检验 法 是 一 种 基于 
样本 频率 进行 统计 检验 的 正 态 检验 方法 ,通常 利用 
正 态 分 布 构造 一 个 gg; 的 分 布 统计 量 (自由 度 n= 
2), 当 设 定 显著 性 水 平 (一 般 取 a=0.05 ) 时 ,分 布 统计 
量 大 于 分 位 数 时 ,接受 假设 H: 总 体 服从 正 态 分 布 ; 
否则 拒绝 假设 。 同 时 ,人 工 增 水 作业 中 一 般 只 研究 
降水 量 的 增加 与 否 ,因此 选择 单 边 1 检 验 。 

当 作 业 期 与 历史 期 的 样本 均 服从 正太 分 布 , 其 
差 值 显 著 性 检验 公式 为 : 

( ( 


e k-1)S} + n-1)S° 1. 1 3 
=O | er) Ea) i 


AP yy TIA EAA A SE EASE YE sk 
n 分 别 为 作业 期 和 历史 期 的 样本 容量 ;5;、5, 分 别 为 
作业 期 和 历史 期 的 样本 标准 差 。 

在 检验 中 ,一 般 取 置信 水 平 a=0.05, 当 由 式 (3) 
计算 出 wi 时 ,表明 人 工 增 水 效果 显著 。 此 时 可 对 
增 水 效果 的 置信 区 间 进 行 估计 ,计算 公式 为 : 

E=y, -yo > (Yr -7,)- bza X 


k- 1)? +(n-1)S? 
K ia } 2 fied 
式 中 :万 为 开展 人 工 增 水 作业 后 样本 的 平均 增加 
HE; yw 为 作业 期 间 假如 不 实施 播 云 作 业 时 的 样本 
自然 平均 值 。 该 式 一 般 取 置信 水 平 (1-a=0.9)。 
由 此 可 得 出 : 


(4) 


R-E -V Yw (5) 
Va Yn 
q=% Xo (6) 
Yok 


式 中 :R 为 开展 人 工 增 水 作业 后 样本 的 相对 增加 率 ; 
d 为 人 工 增 水 效果 检验 功效 ,d 值 越 小 检验 功效 越 
好 ,反之 越 差 。 

当 用 检验 法 对 各 时 期 样本 进行 方差 齐 性 检 
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验 ,得 出 在 显著 性 水 平 c=0.05 时 ,存在 显著 差异 ,此 
时 差 值 的 显著 性 可 用 Welch 检验 法 ,计算 公式 为 : 


z=|y yl/ (SE + Sn (7) 


自由 度 (w) 计 算 公 式 为 : 


vate / 0 
n k(k-1) =- 了 


估计 样本 增值 置信 区 间 的 计算 公式 仍 采用 式 


(k-1)S,+(n-1)S? TT 
(4) ,但 式 中 的 | SS 用 式 (7) 


的 [Sk + Sjn 代替 。 
2 结果 与 分 析 


2.1 统计 变量 样本 的 选取 

为 科学 合理 分 析 巴 音 布 鲁 克 山 区 暖 季 人 工 增 
水 作业 效果 ,未 作业 期 的 样本 数量 应 有 10 a 以 上 。 
其 次 ,人 工 增 水 作业 效果 检验 普遍 使 用 的 统计 指标 
以 降水 量 为 主 ,当选 取 降 水 量 外 的 其 他 样本 时 ,其 
他 样本 应 与 降水 量 相关 显著 。 
2.1.1 降水 量 与 径流 量 的 相关 性 分 析 从 图 2 可 知 ， 
1994—2018 年 暖 季 月 均 降 水 量 的 最 大 值 .最 小 值 以 
及 平均 值 分 别 比 1973—1993 年 增加 了 31.58% 、 
6.71% 和 18.97% ,并 且 降 水 量 与 地 表 径 流量 呈 显 著 
2 HE TE 4H X (R°=0.4743, P<0.001) 。 由 此 说 明 ， 
1994 一 2018 年 是 降水 量 快速 增加 时 期 ,而 径流 量 的 
变化 受降 水 量 影响 明显 ,表明 径流 量 作 为 人 工 增 水 
作业 效果 的 检验 指标 是 可 行 的 。 
2.1.2 模拟 径流 量 与 实测 径流 量 、 降 水 量 的 相关 性 
分 析 五 是 反映 流域 下 垫 面 特征 的 唯一 参数 ,计算 
方式 复杂 ,因此 采用 姚 俊 强 等 3 对 巴 音 布鲁克 山区 


一 一 暧 季 累 积 月 均 降 水 量 
80 (a) 年 际 变化 。 一 历史 期 和 作业 期 累积 月 均 降水 量 均值 


20 
1973 1982 1991 2000 2009 
年 份 


地 表 径 流 分 析出 的 参数 元 =1.21 ,作为 本 文 径流 模型 
中 的 参数 。 

根据 图 3a~c 可 知 ,上 暧 季 月 均 径 流 的 模拟 值 与 
实测 值 均 旦 显著 增长 趋势 (Rw=0.2013,Pw=0.0018; 
Rs=0.1396, Ps=0.0097) , 且 两 者 相关 性 显著 (R= 
0.4758, P<0.001) ,但 实测 值 比 模 拟 值 增 速 大 ,突变 
年 份量 。 从 已 有 研究 可 知 ” 汪 ,20 世 纪 90 年 代 正 是 
西北 干旱 区 冰雪 融 水 量 增加 最 明显 时 期 。 因 此 ,两 
者 增 速 与 突变 年 份 差异 明显 ,与 冰雪 融 水 量 等 其 他 
因素 有 关 。 由 图 3d~f 可 得 ,模拟 径流 量 与 降水 量 的 
波动 性 基本 一 致 ,增长 速率 也 相差 不 大 ,两 者 不 仅 
相关 性 极为 显著 (R=0.9971,P<0.001) ,在 M-K 检 验 
中 变化 趋势 也 几乎 一 致 。 

综 上 可 知 ,选择 实测 径流 量 进行 效果 检验 ,会 
导致 检验 结果 存 有 误差 ,而 基于 Budyko 模 型 估算 出 
的 径流 量 ,不 仅 能 准确 反映 出 降水 量 的 变化 趋势 ， 
还 能 极 好 反映 出 降水 量 对 径流 量 的 影响 。 因 此 , 模 
拟 径流 量 比 实测 径流 量 更 适宜 当 检验 样本 。 

2.2 序列 试验 的 结果 分 析 

由 表 1 可 知 , 暧 季 人 工 增 水 作业 前 后 各 样本 绝对 
增加 的 数值 排序 为 :实测 径流 量 > 降 水 量 > 模 拟 径 流 
量 , 相 对 增加 率 排序 为 :实测 径流 量 (23.28%)> 模 拟 径 
流量 (23.19%)> 降 水 量 (18.98%)。 其 中 ,模拟 径流 量 
比 降水 量 的 绝对 增加 数值 小 1.33, 比 实测 径流 量 的 
相对 增加 率 小 0.09%。 由 此 表明 ,分 析 流 域 径流 量 比 
人 研究 单 站 降水 量 ,更 能 反映 出 人 工 增 水 作业 效果 ,而 
模拟 径流 量 的 检验 功效 不 低 于 实测 径流 量 。 

2.3 不 成 对 秩 和 的 显著 性 分 析 

分 别 将 1973 一 1993 年 21 a 历史 期 和 1994 一 

2018 年 25 a 作 业 期 的 降水 量 .模拟 径流 量 、 实 测 径 
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图 2 1973 一 2018 年 巴 音 布鲁克 山区 暖 季 降水 量变 化 及 其 与 径流 量 的 关系 


Fig. 2 Variation of precipitation and its relationship with runoff in Bayanbulak mountain area in warm seasons form 1973 to 2018 
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(0) 模拟 径流 量 与 实测 径流 量 M-K 检 验 
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图 3 
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注 :Uf .U6、U6 分 别 为 模拟 径流 量 、 实 测 径流 量 、 降 水 量 的 Uf 序列 ;Ub,、Ub,、Ub; 分 别 为 模拟 径流 量 、 实 测 径流 量 、 降 水 量 的 Ub 序列 。 
巴 音 布鲁克 山区 暧 季 模 拟 径 流量 与 实测 径流 量 和 降水 量 的 年 际 变化 ,相关 性 及 M-K 检 验 


Fig.3 Interannual variability, correlation and M-K test of simulated runoff versus measured runoff and precipitation in 


Bayanbulak mountain area in Warm seasons 


表 1 序列 检验 分 析 


Tab.1 Sequence test analysis 


统计 变量 历史 期 均 作业 期 均 绝对 增加 值 
值 (yy,) 值 (yy,) (AE) 
降水 量 /mm 43.53 51.79 8.26 
模拟 径流 量 /ms” 29.88 36.81 6.93 
实测 径流 量 /ms” 152.55 188.07 35.52 


流量 , 按 从 小 到 大 顺序 进行 秩序 列表 计算 ,并 依据 历 
史 期 21 a 的 秩 和 进行 比较 (容量 小 于 作业 期 25 a), 
得 出 暧 季 人 工 增 水 作业 效果 不 成 对 秩 和 检验 表 


4E FA 


(42). HRR FEAR CE LI BY, u (ES) 
于 -1.64, 表 明 巴 音 布鲁克 山区 暖 季 作业 期 的 降水 量 、 
模拟 径流 量 与 实测 径流 量 均 比 历史 期 显著 增加 。 
2.4 t 检 验 的 效果 评估 与 检验 

2.4.1 正 态 分 布 检验 将 各 时 期 降水 量 、 模 拟 径 流 
量 、 实 测 径流 量 分 别 作为 样本 ,利用 S-W 检验 法 ,可 
得 出 不 同样 本 在 各 时 期 的 正 态 分 布 情况 (图 4)。 对 
比 发 现 , 除 作业 期 实测 径流 量 与 历史 期 模拟 径流 量 
的 实际 值 概率 与 理论 概率 较为 离散 外 ,其 余 样 本 均 
较为 相近 ,但 整体 看 各 样本 实测 趋势 与 理论 直线 基 
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表 2 不 成 对 秩 和 检验 分 析 
Tab. 2 Unpaired rank sum test analysis 


历史 期 样本 作业 期 样本 秩 和 统计 检验 显著 性 水 


容量 (mm) ”容量 (ms) (T) (E(u) Æ (P) 
降水 量 21 25 378 -2.56 0.011 
模拟 径流 量 21 25 382 -2.45 0.014 
实测 径流 量 21 25 363 -2.87 0.004 


本 重合 。 运 用 SPSS 27.0 版 本 计算 各 样本 正 态 分 布 
概率 (P) ,得 出 实测 径流 量 Pyw=0.18、Pnw=0.13, 模 
拟 径 流量 Pyw=0.23 .Pw=0.71, 降 水 量 P jw=0.29、 
Ptew=0.87, 均 大 于 显著 性 水 平 a =0.05。 综 上 所 述 ,6 
组 数据 源 服 从 正 态 分 布 。 

2.4.2 方差 检验 进行 方差 齐 性 下 检验 时 ,一 般 常 将 
方差 大 的 作为 第 一 样本 ,方差 小 的 作为 第 二 样本 ， 
通过 比较 1973 一 1993 年 (历史 期 ) 与 1994 一 2018 年 


1.0 F(a) 历史 期 降水 量 。 


降水 量 /mm 


1.0 [ (0) 历史 期 模拟 径流 量 。 “， 
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模拟 径流 量 /m3*s | 


1.0 r (e) 历史 期 实测 径流 量 


(作业 期 ) 各 样本 的 方差 ,得 出 样本 方差 变化 的 检验 
值 ( 表 3)。 给 定 信和 度 a =0.05 WF, HW Foos EK 
于 降水 量 与 模拟 径流 量 的 下 值 ,小 于 实测 径流 量 的 
下 值 ,表明 除 实测 径流 量 的 总 体 方差 差异 性 显著 外 ， 
其 余 样 本 的 总 体 方差 未 改变 。 因 此 ,降水 量 与 模拟 
径流 量 适合 :检验 法 ,而 实测 径流 量 采 用 Welch 检 
验 法 。 

2.4.3 检验 计算 由 表 4 可 知 , 降 水 量 的 检验 功效 最 
好 ,作业 后 只 需 增 值 为 自然 均值 的 11.59% ,就 能 ! 
著 检测 出 效果 ;模拟 径流 量 的 检验 功效 次 之 , 相 比 
实测 径流 量 的 检验 功效 值 低 3.72%。 在 显著 性 水 平 
90% 的 置信 区 间 ,模拟 径流 量 的 相对 增 率 比 实测 径 
流量 少 3.78% , KEKEE 2.65% ,考虑 到 实测 径流 
量 受 冰雪 融 水 等 影响 ,以 及 整个 汇 水 区 降水 与 单 站 
降水 的 差异 性 ,模拟 径流 量 的 检验 结果 可 信和 度 更 
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图 4 降水 量 模拟 径流 量 和 实测 径流 量 在 历史 期 与 作业 期 的 正 态 分 布 


Fig.4 Normal distributions of precipitation, simulated runoff and measured runoff in historical period and operation period 
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表 3 方差 齐 性 下 检验 
Tab.3 F-test for homogeneity of variance 
时 其 降水 量 模拟 径流 量 实测 径流 量 自由 度 Fons 
方差 F 方差 F 方差 F n, n, ( œ =0.05) 
1973—1993 年 38.680 1.997 28.615 1.889 382.805 4.324 20 24 3.226 
1994—2018 4E 77.241 1.997 54.051 1.889 1655.103 4.324 20 24 3.226 
表 4 zt 检 验 分 析 
Tab.4 t-test analysis 
统计 变量 ”检验 功效 (d)/% t z ”自由 度 (v) i(a=0.1) ws(a=0.2) 显著 性 水 平 (P) 相对 增加 量 (Z) 相对 增 率 (R)/% 
降水 量 11.59 3.493 - 45 1.301 0.681 0.0006 5.38 mm 12.05 
模拟 径流 量 14.04 3.456 一 45 1.301 0.681 0.0006 4.53 ms 14.70 
实测 径流 量 17.76 一 3.985 93 1.291 0.677 0.0001 28.48 ms 18.48 


好 。 综 上 可 得 ,基于 Budyko 模 型 得 出 的 径流 量 不 但 
检验 精度 较 高 , 且 有 助 于 提升 检验 功效 。 同 时 , 巴 
音 布鲁克 山区 开展 人 工 增 水 作业 后 ,降水 量 与 径流 
量 增 加 效果 显著 。 


3 讨论 


由 于 巴 音 布 鲁 克 山 区 长 时 序 降水 数据 多 为 单 
站 观测 数据 ,而 人 工 增 水 影响 区 域 又 较为 广阔 ,只 
利用 单 站 降水 数据 开展 增 水 效果 分 析 会 一 定 程 度 
上 存在 影响 。 因 此 ,基于 地 表 径 流 由 整个 汇 水 区 降 
水 形成 理论 ?22 ,以 及 作业 区 地 表 径 流 与 代表 性 站 
点 降水 相关 性 极为 显著 (P<0.001) 的 结果 ,选择 径流 
量 作为 增 水 效果 又 一 检验 指标 ,具有 一 定 合 理性 。 
但 山区 径流 形成 复杂 ,直接 用 于 效果 检验 将 导致 结 
果 偏 差 ,而 Budyko 模型 能 较 好 地 解决 该 问题 ,如实 
测 径 流量 的 相对 增 率 比 降水 量 多 达 6.43% ,可 模拟 
径流 量 不 仅 相 对 增 率 与 降水 量 仅 差 2.65% ,检验 功 
效 值 也 比 实测 径流 量 低 3.72%。 由 此 说 明 , 除 降 水 
作为 检验 指标 外 ,基于 Budyko 模 型 的 径流 量 作为 检 
验 指标 也 具有 较 好 的 精度 与 检验 功效 。 

人 工 增 水 实际 作业 中 , 受 安全 、 人 为 等 因素 影 
响 ,作业 区 往往 位 于 偏远 、 少 有 探测 设备 的 区 域 。 
因此 ,在 常年 作业 区 少 有 降水 观测 站 点 时 ,利用 与 
降水 显著 相关 的 其 他 统计 量 作为 检验 指标 ,有 助 于 
初步 回答 作业 效果 问题 。 在 本 研究 中 ,单个 降水 测 
站 实际 只 能 观测 部 分 作业 区 域 , 而 Budyko 模 型 受气 
候 弹性 系数 唯一 参数 五 影响 , 易 导 致 模拟 径流 量 存 
在 差异 ”。 因 此 , 巴 音 布鲁克 山区 人 工 增 水 作业 效 


注 :t\z 为 统计 变量 计算 出 的 显著 性 检验 值 ;hi ;为 增 水 效果 检验 置信 水 平 值 。 


果 应 在 降水 检验 结果 之 上 ,在 模拟 径流 量 检验 结 
之 下 , 即 相 对 增 水 率 在 12.05%~14.7% 之 间 。 此 外 ， 
数理 统计 检验 结果 受 方 法 ,参数 的 影响 ,一 般 只 具 
参考 性 ,而 从 物理 机 制 变化 开展 的 物理 检验 结 

更 具 说 服 性 。 后 续 研 究 中 ,将 利用 天 气 雷 达 回 波 等 
言 息 分 析 山 区 人 工 增 水 作业 效果 ,以 期 验证 统计 检 
验 的 结果 。 
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4 结论 


(1) 由 于 巴 音 布鲁克 山区 作业 区 范围 广阔 ,长 
时 序 代表 性 降水 观测 站 点 数量 较 少 , 想 要 整体 检 
验 出 人 工 增 水 作业 效果 ,借助 与 降水 显著 相关 的 
开 都 河 地 表 径 流量 作为 统计 检验 指标 (R=0.4743， 
P<0.001) ,是 具有 一 定 可 行 性 的 。 

(2) 地 表 径 流 构 成 复杂 ,直接 用 于 人 工 增 水 效 
检验 易 导 致 结果 偏差 ,通过 Budyko 模型 得 出 的 
径流 量 ,不 仅 与 实测 径流 量 的 变化 趋势 具有 同步 
性 ,还 与 降水 量 的 相关 程度 接近 于 1。 因 此 ,基于 
Budyko 模 型 的 径流 量 主要 代表 降水 对 地 表 径 流 的 
影响 量 , 相 比 实测 径流 量 更 适合 作为 检验 指标 。 

(3) 在 序列 试验 与 1 检验 中 ,模拟 径流 量 与 实测 
径流 量 的 相对 增 率 均 大 于 降水 量 ;模拟 径流 量 较 实 
测 径 流量 的 检验 功效 值 低 3.72%。 表 明 整 个 汇 水 区 
降水 与 单 站 降水 具有 差异 性 ,而 选择 Budyko 模型 的 
径流 量 作为 检验 指标 ,不 仅 检验 精度 有 提升 ,结论 
也 较为 合理 。 

(4) 选取 可 信和 度 为 90% 的 置信 水 平 , 得 出 暖 季 
月 均 实测 径流 量 增加 了 28.48 més ,相对 增 率 为 


Nm 


1970 


FHEWY 


46% 


18.48%; 降 水 量 增加 了 5.38 mm ,相对 增 率 为 12.05% ; 


模拟 径流 量 增加 了 4.53 mes 


,相对 增 率 为 14.7% 。 


证 明 巴 音 布鲁克 山区 暧 季 开 展 人 工 增 水 作业 后 , 降 


水 量 及 径流 量 


EE. 


的 增加 效果 显著 。 
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Effect of artificial precipitation enhancement in Bayanbulak mountain 


area in warm seasons based on Budyko model 


DIAO Peng, LI Gang, YUAN Xianlei, WEN Chun 


(Meteorological Bureau of Bayingol Mongolian Autonomous Prefectrue, Korla 841000, Xinjiang, China) 


Abstract: Although precipitation serves as a common statistical indicator, its utility is limited by the availability 
of representative station data, which is affected by geographical, economic, technical, and other effects. Such limi- 
tations can introduce inaccuracies in the assessment of the effectiveness of artificial precipitation enhancement. 
To address these challenges, this study focused on daily meteorological data from the Bayanbulak meteorological 
station and monthly runoff data from the Dashankou hydrological station in Xinjiang, China, from May to Sep- 
tember in 1973—2018. This study involved building a runoff simulation equation using the Budyko model. To ob- 
jectively and quantitatively analyze the impact of artificial precipitation enhancement in the Bayanbulak moun- 
tain area during warm seasons, sequence tests, unpaired rank sum tests, and t-tests were used. The results show 
the following: (1) The use of the Budyko model for building simulated runoff not only synchronized with the 
changing trends and growth rates of precipitation but also demonstrated a highly significant correlation (R= 
0.9971, P<0.001). This shows that the simulated runoff not only captured the overall precipitation trends but also 
quantified the impact of precipitation on runoff. (2) Utilizing an unpaired rank sum test and t-test, it was found 
that both precipitation and runoff significantly increased (P<0.02) after artificial precipitation enhancement when 
considering measured runoff, simulated runoff, and precipitation as statistical variables. (3) The best statistical in- 
dicator for assessing the impact of artificial precipitation enhancement was precipitation, accounting for only 
11.59% of the added value to statistically test the effect. Compared with the measured runoff, the test efficiency 
value of the simulated runoff decreased by 3.72%, indicating that the test efficiency was improved. (4) With a se- 
lected significance level of a 90% confidence interval, the absolute increase in precipitation was 5.38 mm, repre- 


' indicat- 


senting a relative increase rate of 12.05%. The absolute increase in the simulated runoff was 4.53 m°» s` 
ing a relative increase rate of 14.70%. The absolute increase in the measured runoff was 28.48 m*+s ', correspond- 
ing to a relative increase rate of 18.48%. This indicates the important impact of artificial precipitation enhancement 
during the warm seasons in the operation period (1994—2018) compared with the historical period (1973—1993) in 
the Bayanbulak mountain area. 


Key words: artificial precipitation enhancement; statistical analysis; effect evaluation; Bayanbuluk mountain area 


